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Sunyaev-Zeldovich signal processing and temperature-velocity 
moment method for individual clusters 
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ABSTRACT 

Future high resolution, high sensitivity Sunyaev-Zeldovich (SZ) observations of individual 
clusters will provide an exciting opportunity to answer specific questions about the dynam- 
ical state of the intra-cluster medium (ICM). In this paper we develop a new method that 
clearly shows the connection of the SZ signal with the underlying cluster model. We include 
relativistic temperature and kinematic corrections in the single-scattering approximation, al- 
lowing studies of hot clusters. In our approach, particular moments of the temperature and 
velocity field along the line-of-sight determine the precise spectral shape and morphology of 
the SZ signal. We illustrate how to apply our method to different cluster models, highlight- 
ing parameter degeneracies and instrumental effects that are important for interpreting future 
high-resolution SZ data. Our analysis shows that line-of-sight temperature variations can in- 
troduce significant biases in the derived SZ temperature and peculiar velocity. We furthermore 
discuss how the position of the SZ null is affected by the cluster's temperature and velocity 
structure. Our computations indicate that the SZ signal around the null alone is rather insen- 
sitive to different cluster models and that high frequency channels add a large leverage in this 
respect. We also apply our method to recent high sensitivity SZ data of the Bullet cluster, 
showing how the results can be linked to line-of-sight variations in the electron temperature. 
The tools developed here as part of SZpack should be useful for analyzing high-resolution SZ 
data and computing SZ maps from simulated clusters. 
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1 INTRODUCTION 



Since the e arly measuremen t s of t h e SZ e ffect in the 80's and 
IBirkinshaw et all Il984 1 19911 : lLamarre et all 1 19981 : 



90's (e.g., 

iHughes & Birkinshawlll998l) the observational possibilities in the 
microwave band have evol ved at an impressive rate. T oday the 
thermal SZ (thSZ) effect dZeldovich & Sunvaevi Il969l) is rou- 
tinely detected for several hundred galaxy clusters out to redshift 
z — 1 ^Benson et alj[20 04: Mar riage et alj|2oTll : Iwilliamson et alj 
1 201 it IPlanck Collaboration et al.1 1201 lh . and in some cases evi- 
dence for_Jhe_rjresence_of the s maller kine matic SZ (kSZ) ef- 
fect (Sunyaev & Zeldovich 1980]) is found I Benson et al J 20031: 



iKorngut et alJl201ll : lMroczkowski et alj|2012l : lHand et alj|2012h . 

SZ cluster observations are a pow e rful cosmological tool 
(see iRephaelil Il995d : IBirkinshaw! 1 19991 : ICarlstrom et ail 120021 : 
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iKomatsu & SeliatJ 12002! , for overview), but to realize their full 
constraining power for precision cosmology it is critical to un- 
derstand the structure and evolution of the intra-cluster medium 
(ICM) in more detail. Already now a number of high-resolution 
SZ experiments, including ALM41 CARM^E CCAlfl and 
MUSTANC0, are underway or planned, promising a dramatic 
increase in sensitivities, spatial resolution, and spectral coverage 
over the next few years. Current generation high-resolution SZ 
observations of individual clusters have already revealed rich 
phenomena (e.g., shocks, substructures, re lativistic partic l es) in 
the atmospheres of merging clusters (e.g., 
iKitavama et alj 120041 : [Colafrancesc o et alj 



Komatsu et al]|200ll: 



2011: Korngut et all 



1 Atacama Large Millimeter/submillimeter Array 

2 Combined Array for Research in Millimeter-wave Astronomy 

3 Cornell Caltech Atacama Telescope 

4 Multiplexed Squid TES Array at Ninety GHz 
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Menanteau et all |2012|; iMroczkowski et al] |2012| ; 
Zemcov et alj2012l ; |Prokhorov & Colafrancescoll2012l) . Upcoming 



SZ experiments should enable a host of exciting measurements 
of important physical processes that shape the properties of the 
ICM and their evo l ution, including the electr on temperature 
jPointecouteau et all 1 19981 : lHansen et alj l2002r) . the peculiar 
veloc ity of the cluster, internal bulk and turbulent gas motions 
(e.g. IChluba & Mannheiir]|2002l ; iNagai et alj|2003l ; ISunvaev et all 
l2003| - |DiegoetalJl2003h . and non- equilibrium electrons produced 
by merger and accr etion shocks dMarkevitch & Vikhlininl |2007| ; 
iRudd & NagaHl2009h . 

These encouraging prospects also raise a number of impor- 
tant problems that must be addressed before the rich information 
contained in the future SZ data can be fully exploited. One is sim- 
ply related to the precise and fast computation of the SZ signal 
given basic parameters of the scattering medium, such as the Thom- 
son scattering optical depth, r, the electron temperature, T t , and 
bulk velocity, j8 c , while accounting for relativistic temperature and 
kinematic correctio ns. Previously, this issue has been addressed 
by several groups (IRephaelil 1 1 995 bl; [Challinor & Lasenbvl Il998 l; 
Itohetalj[l998; Sazonov & SunvaevHl998l; iNozawa et al.lll998tj; 



Challinor & Lasenbvl ll999; Chluba et al. 2005?) by means of Tay- 



lor expansions for the SZ signal. While the evaluation of these 
expansions is very fast, they are limited to rather low tempera- 
ture gas (Fig. [2] shows that for T c > 13keV this approach breaks 
down). One alternative is direct numerical integration of the Boltz- 
mann collisi on term making use of the symmetries of the scatter- 
: problem JWrighdl979l;lDolgov et al.l200ll;lNo"zawa et al.l2009l ; 



ing proble m HWrigtigiy /ytlDotgov et al.l2UUll ; IJNozawa et al.lZUU^ ; 
IPoutanen & VurmfeOlOh ." but this is time-consuming and not well- 



suited for extensive parameter estimations or computations of the 
SZ signal from simulated clusters. On the other hand, a fast but not 
as precise and flexible approach is simple tab ulation of the SZ sig- 
nal or the use of convenient fitting function 1 Nozawa et al JI2000I ; 
Itoh&Nozawa 2004; Shimon & RephaebfeOQ^ . 

Recently, Ichluba et alj J2012l. CNSN in the following) devel- 
oped a method in the middle of these extremes. In their work, a new 
set of frequency-dependent basis functions was computed numeri- 
cally to allow very fast and precise calculation of the SZ signal. The 
basis functions are informed by the underlying physics of the scat- 
tering problem and thus are ideally suited for future SZ signal pro- 
cessing and parameter estimation. The associated routines are part 
of SZpAClfl However, several extensions are required. Firstly, so 
far high precision (relative accuracy =2 0.001%) calculations with 
SZpack were limited to T c < 25 keV. This problem can be easily 
overcome using the method of CNSN by appropriate extensions of 
the basis functions, as we explain in Sect. 12.41 With SZpack vl.l, 
which is presented here, it is now possible to compute the SZ signal 
for T e < 75 keV and /3 C < 0.01 to =; 0.001% relative precision at 
practically no computational cost. This precision and range of pa- 
rameters covers all physically relevant cases and hence provides an 
important preparation for SZ parameter estimation without signifi- 
cant limitations. 

Secondly, line-of-sight variations of the temperature and ve- 
locity field (with any of the aforementioned methods) can only 
be accounted for by means of additional 1 -dimensional integrals; 
this again makes extensive SZ parameter estimation expensive. Es- 
pecially when computing the SZ signal from cluster simulations 
the problem becomes very demanding, even if evaluation for single 
gas parameters (r, T e , and /? c ) is extremely fast. Here we reformu- 
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late the representation of the SZ signal to overcome this limitation. 
We utilize that in the single-scattering approximation, frequency- 
dependent terms can be separated from temperature- and velocity- 
dependent contributions (Sect. [2)- This implies that the SZ signal 
for a given cluster model can be calculated using appropriate mo- 
ments of the temperature and velocity field. While this means that 
a finite number of 1 -dimensional integrals along different lines-of- 
sight have to be evaluated, this separation still greatly reduces the 
computational burden because afterward the SZ signal at any fre- 
quency can be computed as simple matrix multiplication. 

While the method developed in Sect. [2] is both precise and 
fast, delivering quasi-exact results for the SZ signal through differ- 
ent lines-of-sight for any cluster atmosphere, in the future analysis 
of high resolution, high sensitivity SZ data another simplification 
is possible. In Sect. [3] we show that for typical cluster models the 
smoothness of the temperature and velocity profiles allows mini- 
mizing the number of parameters needed to accurately describe the 
SZ signal, resulting in a second set of moments that are related to 
the line-of-sight temperature and velocity dispersions and higher 
order statistics (see Eq. ([18} for instance). The associated expan- 
sion of the SZ signal around the mean becomes perturbative and 
the number of moments needed to describe the SZ signal depends 
directly on the observational sensitivity. 

Our formulation furthermore allows direct separation of 
frequency-dependent from spatially varying terms, providing a 
clear link between morphological changes of the SZ signal and 
cluster parameters. For example, the presence of large-scale, 
post-merger cluster rotation can introduce a bipolar kSZ signal, 
which is related to a spatially varying average line-of-sight ve- 
locity. The superposition of thSZ with this rotational kSZ implies 
small frequency-dependent changes of the clusters morphology 
dChluba & Mannheinj [20021) . Similarly, variations of the electron 
temperature along the line-of-sight intr oduce morphological e ffects 
(as also pointed out more recently by IProkhorov et al.ll201 lh . and 
as we explain here, spatial variations of temperature and velocity 
moments are the source of these morphological changes. The mo- 
ments therefore constitute the main observables of high-resolution, 
high-sensitivity SZ observations and their interpretation is the main 
challenge for future SZ parameter estimation and in the reconstruc- 
tion of the cluster's temperature and velocity structure. 

Armed with these tools, we address a number of questions 
that are related to the effect of temperature and velocity variations 
on the SZ signal, with particular focus on parameter degeneracies 
and instrumental aspects. For example, we explicitly discuss the ef- 
fect of angular resolution and frequency filters on the SZ signal, as 
well as different corrections to the location of the SZ null. All these 
aspects, if ignored, lead to biases in the deduced cluster parame- 
ters. We illustrate this for several examples, using both mock SZ 
measurements as well as recent SZ data. We furthermore develop 
several tools for SZ parameter estimation which are now part of 
SZpack. These should be useful for computing the SZ signal from 
cluster simulations and in the analysis of future high resolution, 
high sensitivity SZ measurements. 



2 COMPUTING THE SZ SIGNAL USING 

TEMPERATURE AND VELOCITY MOMENTS 

In this section we introduce the new temperature-velocity moment 
method to compute the SZ signal for general cluster atmospheres. 
This section is rather technical and mainly for readers interested in 
the computational details. 
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The SZ effect is caused by the scattering of CMB photons by 
moving electrons. For a small volume element of scattering elec- 
trons the SZ signal only depends on the electron temperature, T e , 
their total bulk velocity, j8 c , the direction cosine of this velocity with 
respect to the line-of-sigh{| /i c = jS c • y, and the Thomson optical 
depth of the scattering volume element, At. Both T e and At are de- 
fined in the rest frame of the scattering volume element, while j8 c is 
defined with respect to the CMB rest frame. 

For now, we shall assume that the observer is at rest in the 
CMB frame. The conversion of the aforementioned parameters to 
the SZ signal can be expressed as A/(x) = AtFsz(x, T t ,p c ,ji c ), 
where F sz is a non-lin ear function and x = hv/kT with CMB tem - 
perature T = 2.726 K fFixsen et alJl996l . lFixsen & Mafheill 20021) . 
The change in the CMB intensity can be further rewritten as 

A/(x) « At /„ x 3 [Y(x, T e ) + p\M(x, T c ) 

+/? C P, 0u c ) C(x, r.) + fiPiiflc) Q{x, T e )\ , (1) 

where I = (2h/c 2 )(kT /h) 3 « 270 MJy sr" 1 and P,(x) denote Leg- 
endre polynomials. The term Y(x, T e ) describes the purely thermal 
SZ effect with temperature corrections included, while terms oc /3 C 
are related to kinematic effects. 

From previous analysis of the SZ effect, it is furthermore clear 
that the functions Y(x, T e ), M(x, T e ), D(x, T e ) and Q(x, T e ) can all 
be described using an appropriate set of frequency-dependent basis 
functions and temperature-dependent coefficients, where the latter 
encode the spatial dependence. The integrated SZ signal along a 
given line-of-sight is therefore determined by appropriate moments 
of the cluster's temperature and velocity field. These contain the 
desired information about the cluster gas and structure; as such they 
define the observables of the SZ measurement, and the aim will be 
to use the moments to learn about the cluster gas. 

2.1 Low temperature gas (kT e < lOkeV) 

As discussed in IChluba et alj <2012h . for electron gas tempera- 
tures kT e < 5keV - lOkeV an a symptotic expansion of the 
Boltzmann colli sion term (similar to jchallinor & Lasenbylll99£ ; 



Nozawa et al 



1998a 



Itoh et alj|l998l; ISazonov & Sunvaevlll998l : _ 
Challinor & Lasenbvlll99^ can be used to represent the SZ signal 
with high precision. From Eq. (25) of CNSN, up to some speci- 
fied correction order, fe max , in the electron temperature, it therefore 
follows (see Appendix [A] for more details): 



r low (x,r e ) = JV,(x)< 



M 



"JUilA 

v (x,r e ) = -AfW + ^Mfwf 



D , ^(x,T c ) = g(x)+YD 1 r(x)e k e + 



{fix, r e ) = H Q(x) + £ el ow (x) e? 1 



(2a) 



(2b) 



(2c) 



(2d) 



k=0 

with 8 C = kTJm e c 2 , M(x) = Y (x) + @{x\ Q(x) = xe*/[e* - if, 
and Q = x{jcoth(x/2). The functions Yt are defined as in CNSN, 
while M' cw , D l ™, and g[ ow are given by Eq. {A3). To give an exam- 
ple, Y = Q - AQ describes the usual (non-relativistic) thSZ effect 

6 In the following bold font denotes 3-dimensional vectors and an addi- 
tional hat means that it is normalized to unity. 



dZeldovich & Sunvaevl 1 19691). while the term oc Q is related to the 
kSZ effect dSunvaev & Zeldovicd[l98(ih . Both M(x) and Q(x) de- 
scribe the lowest order kinematic terms oc fi 2 . 

For a fixed line-of-sight through the cluster medium, the to- 
tal SZ signal is determined by integration over At. With the de- 
composition given above it is convenient to introduce the following 
line-of-sight temperature and velocity moments: 



/) = dT, 



b m = 



r.dr, 



bf = J ^P 2 0i c )^dr, 



(3) 



where the integrals are carried out in the cluster frame, with the 
condition T c < r e ,iow for k > and for y m . Here y {k) denotes 
the generalized y-parameters, while bf } take into account the ef- 
fect of the clusters global and internal gas motion. For example, 
y (0} = j 8 e dT is the usual line-of-sight y-parameter or average ther- 
mal pressure of the electrons, while bf* (no temperature depen- 
dence) is proportional to the average velocity of the cluster medium 
along the line-of-sight weighted by the electron number density. 
The optimal values for r c j ow and k mlix depend on the required pre- 
cision and will be specified below (see Sect. l2.4t . With Eq. ([3) we 
can now define the moment vector 



y 

bj 



(y T ,b],b T ,b T 2 ) 

(/>,..., y*—)) 
{bf\ 



,b ikm 



(4) 



where we arranged the entries of wii ow with respect to orders in 
/J c . Notice that the dimensions of the vectors bj in principle can 
differ from k max + 1. In particular, for the velocity moments usually 
fewer terms in the electron temperature are required to describe the 
SZ signa l accurately, since th ey only lead to very small corrections 
(e.g., see lChluba et alj2012h . 

If we assume that the SZ signal is observed at m frequencies, 
{x,}, we can furthermore introduce the signal vector for the contri- 
bution of the low temperature gas, S^ w = (A/ tow (xi), A/i ow (x„,)). 
This defines the matrix equation 

*^low ~ Fl ow /Wl ow , (5) 

where each row of the matrix F) ow reads 

(Flow), = /„ 4 (fL. 9, DL, \m, Mj ow , Hq, gfj , (6 ) 



with Y' low = (F 



n„,,<), D' low = (£#«, and so on. For 

the i' h row of Fi ow the frequency-dependent functions evidently 
have to be evaluated at the required x,. The SZ signal caused by 
low temperature gas can therefore be computed as a simple matrix 
operation once the low temperature- velocity moment vector, m\ ow , 
is known. The columns of the moment matrix, F low , are the values 
of the basis functions at the required frequencies. 



2.2 High temperature gas ikT e > lOkeV) 

For gas with temperatures kT c > lOkeV the convergence of the 
asymptotic expansion given above becomes slow (see Fig. |2). 
However, recently CNSN found an alternative set of frequency- 
dependent basis functions that allow very accurate description of 
the SZ signal up to high temperatures and bulk velocities. The con- 
vergence radius of the CNSN expansion around the chosen pivot 
temperature (in CNSN kT e fj = 15keV was used) is also limited 
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(see Fig. |2), but in combination with the low temperature expan- 
sion it allows covering a large part of parameter space. 

Like for the asymptotic expansion, the signal is determined 
by particular temperature and velocity moments, but this time the 
weighting differs slightly from those of Eq. ((3}. With the expres- 
sions given in CNSN it is straightforward to show (see Appendix [B] 
for more details) that up to some specified order, fc,^ ax , of the elec- 
tron temperature one has 



M h '"\x, T e ) = £ Mf E \x) N(6 e ) 9 k 

k=0 
Km 

D tu " h (x,T e ) = Y J D l l ,gh ( X )N(e c )e k e 

k=0 
Km 



(7a) 



(7b) 



(7c) 



(7d) 



[1 " H +^ + 0(flg)]. Here, 



withJV(0 c ) = t^t^ - ^ 

K 2 (x) denotes the modified Bessel functions of second kind. The 
functions Yf gh , M*' gh , d" 8 ", and g* igh are defined by Eq. (|B3j. All 
temperature-independent kinematic terms were already taken into 
account by Eq. l(2j, so that they do not reappear here. Also, in gen- 
eral k* mxi t & max , although below we usually set k^ ax = 

In analogy to the low temperature gas case we introduce the 
following line-of-sight temperature and velocity moments: 



N(9 c )ffL&r. 



I 

f APiOieWCftO^dr, cf = J #:P 2 (^)JV(A)^dT, 



where the integrals are carried out in the cluster frame, with the 
condition r e j ow < T c < r e _hi g h- The optimal value for r e _hi g h depends 
on the required precision. Also, one can split the temperature range 
up into different parts, each with their own set of basis functions 
defined on the intervals T c e (r^! gh , r[ high ], as will be specified in 
Sect. 12.41 With this we define the moment vector 



Vgh = (z ,ci,e ,e 2 ) 
z r = (/»,..., z«->) 



(9) 



cf = (<f. 



The signal vector for the contribution of the high temperature gas, 
Shi g h> is defined similar to Si ow , and it can be obtained with 

Shigh = Fhigh"*high> (10) 

where each row of the matrix Fhi a h reads 



(Fhigh)i - h x \ (^high' ^\igh> ^high> fihigh) i 



(ID 



with F h r lgh 



high 



/high 



gh 



(D„ gh , 



D'°, w ), and so 



on, all as above. Again the SZ signal caused by high temperature 
gas can be expressed as a simple matrix multiplication, once the 
moments are determined. This reduces the computational burden 
to calculation of the temperature-velocity moment vector which 
solely depends on the cluster atmosphere. 

2.3 Total line-of-sight SZ signal 

With the definitions of the previous sections, the total SZ signal is 
given by S = Si ow + Shi e h- Introducing the total cluster temperature- 



velocity moment vector, m T = (pt[ ov/ , m l igh )< an d the frequency- 
dependent matrix, F = (Fi ow , F Mgh ), one has 

S = Fm. (12) 

It is clear that the dimension of S defines the maximal number of 
moments that could possibly be deduced from the SZ data. How- 
ever, in the presence of noise, foregrounds, and correlations be- 
tween the moments, one naturally has dim(wi) < dim(S). 

It is furthermore important that for a given experimental 
precision, an optimal combination of the basis functions can be 
found which minimizes the number of moments required to ac- 
curately represent the SZ signal. In particular, with the approach 
of CNSN one can vary the reference/pivot temperature, r e0 , and 
number of reference points used in the computation of the basis 
functions to improve the temperature coverage of the approxima- 
tion. In that case the moment vector can be cast into the form 



{m. 



T T T 

low 1 high!' highH' 



.) with different temperature regions, 



[0, r c , low ], (7V, ow , 7^], (7^ hlgh , r° Ugh L (r° Ugh , rg^J, and so on. 
Here the subscript 'low' is used to indicate that the asymptotic ex- 
pansion is applied for those moments, while for any moments with 
subscript 'high' we formulate the basis using CNSN. We will dis- 
cuss the associated optimization problem in Sect. 12.41 

With the formulation given above it is also straightforward to 
include the effect of angular resolution and frequency bands on the 
SZ signal. The effect of angular resolution is accounted for by spa- 
tially averaging the temperature-velocity moments, i.e., m — » (m), 
where (...) denotes angular/spatial average. The bandpass can be 
taken into account by means of a matrix W. With this the SZ signal 
in more general can be expressed as 



S = W F (m) + n, 



(13) 



where we also added noise to the problem. In a similar way possible 
contaminations by (spatially) smooth foregrounds, radio sources, or 
dusty-star-forming galaxies (DSFGs) can be incorporated. 

One of the benefits of the moment method described here is 
that for a given set of frequencies the moment matrix only has to be 
computed once. This, for example, makes computation of the SZ 
signal from simulated cluster very efficient and accurate. However, 
Eq. dl2l > is still mainly interesting from the computational point of 
view because the entries of the moment vector are not independent. 
For instance, all moments related to are non-negative and one 
also expects y {k> > y {k+i) . This imposes rather complicated priors 
and correlations among the different entries of the moment vector 
with the actual dimensionality of the problem being much smaller. 
We will show below that for future SZ observations only a few pa- 
rameters are required to accurately describe the SZ signal, although 
the number of entries in the moment vector, m, is much larger. 

Finally, we mention th at the effect of the obser vers motion 
with respect to the cluster dChluba etal.ll2005ll2012l) can be in- 
cluded using simple Lorentz-transformation of the frequencies and 
corresponding angles into the CMB rest frame, to account for the 
effect of Dop pler boosting and relativistic light aberration (e.g., 
IChlubalfeoill) . However, for the discussion below these aspects of 
the problem are not crucial and will be omitted. 



2.4 Minimizing the required number of moments 

One of the important questions is how many moments are needed to 
describe the SZ signal accurately for a given experimental sensitiv- 
ity and range of gas temperatures. Here it is particularly interesting 
to try minimizing the total number of moments that are required to 
achieve an optimal representation of the SZ signal. To answer this 
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Figure 1. Deviation of the approximation from the numerical result. The 
asymptotic expansion is used to represent the SZ signal. For each of the 
curves, k + 1 temperature orders were included. The departure is expressed 
in units of A/ s = 0.013 MJy sr" 1 and t = 0.01 was assumed. 












k even \ \ \ \ \ ' \ i 

k odd \ * \ 1 \ > \ ' 





T [keV] 

Figure 2. Range of convergence for different approximations. The upper 
panel shows the results obtained with the asymptotic expansion, while for 
the lower panel the basis CNSN with 8 e q = 0.03 was used. The maximal 
departure for 0.1 < x < 30 is expressed in units of A/ s ^ 0.013 MJy sr~' 
and t = 0.01 was assumed for all cases. 



question we first define a fiducial sensitivity for comparison. We 
shall use the kSZ signal of a cluster with line-of-sight optical depth 
r = 1(T 2 and /3 cJ u c = 1(T 3 close to the thSZ crossover frequency 
y c = 217 GHz (x c = 3.83) as benchmark; this gives a distortion with 
amplitude A/ s a 9.76T/? cy u c (fcT ) 3 /[/? 2 c 2 ] sr 1 0.013 MJy sr 1 . 
For an isothermal cluster with kT c m 5 keV electrons and t = 10~ 2 
the maximal thSZ signal is roughly A/ th = 0. 17 MJy sr -1 at x ^ 6.7. 
Therefore, A/ s corresponds to = 8% precision on A/ th . 

Before carrying out additional computations we extended the 
basis of CNSN with additional pivot points (at 6 e = 0.01 and 0.1) 
such that the SZ signal can be represented in a wider range of tem- 
peratures. We provide this basis both in the cluster rest frame and 
the CMB frame. With the current version of SZpack a = 0.001% 
precision is achieved at frequencies 0.01 < x < 30, temperatures 
kT e < 75 keV, and for f} c < 0.01 basically at no additional com- 
putational cost. We furthermore included the necessary database 
directly into SZpack such that no time is consumed loading data. 
In the current implementation evaluation of the SZ signal at 400 
frequencies takes about 0.01 seconds on a standard laptop. These 
routines can also be directly invoked from Python. 

One can now calculate how accurately the different sets of ba- 
sis functions describe the SZ signal for varying k ma . For the asymp- 
totic expansion we show two examples in Fig. [Tj At low tempera- 



tures (upper panel) only a few terms in the expansion are needed to 
achieve a very accurate representation of the SZ signal. The mis- 
match is usually largest at high frequencies, x - 10, while below 
the crossover frequency higher order temperature terms are small, 
even for larger electron temperatures. The lower panel of Fig.Q]in- 
dicates that at higher temperatures the convergence of the asymp- 
totic expansion becomes sl ower, a problem that is well-known from 
previous analysis (e.g., see lltoh et al.ll 1998h . 

A simple calculation can be used to further quantify the con- 
vergence rate of the different basis functions: for a given order in 
temperature we compute the maximal deviation of the approxima- 
tion from the numerical result in the frequency range 0.1 < x < 30. 
For both the asymptotic expansion and the basis CNSN with refer- 
ence temperature G e = 0.03 the results are shown in Fig. [2] One 
can see that for the asymptotic expansion the agreement with the 
numerical result does not improve above kT c a 13keV. At higher 
temperatures the expansion of CNSN performs much better. In par- 
ticular, the convergence radius increases strongly when including 
the first few temperature corrections. One can also observe that for 
the CNSN basis functions with pivot temperature 9 s = 0.03 con- 
vergence above kT c ^ 40keV is not achieved. However, this can 
be overcome by adding another set of basis functions with pivot 
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Table 1. Optimal distribution of temperature pivots, T' a , and regions, 



[0, ^ciowl (r c ,iow, r' high ] 



(^ehi h'^ehiglJ' etc ' ^ or gi ven accuracy goal, 



A//A/ s . At temperature T c < 7* e i ow the asymptotic expansion is used, while 
above expressions based on CNSN are applied. As fiducial accuracy value 
we used A/ s =s 0.013 MJy sr~' and optical depth r = 0.01. Numbers marked 
with asterisk are only estimated upper bounds, although the approximation 
is much better. The number of temperature correction terms in each region 
is itmas. while k tot gives the total number of temperature terms for all re- 
gions. The required number of temperature-velocity moments depends on 
the settings for kinematic corrections and is not further specified here. 



A//A/ s 


^max / ^tot 


[keV] 


e.high e.O 

[keV] 


e,high ' e,0 

[keV] 


j-iii /j-iii 

e.high e,0 

[keV] 


1 
1 
1 


2/12 
3/12 
4/10 


9.1 
12.5 
14.3 


21.6/14 
48.5/25 
75/35 


42/30 
90* /80 


75/55 


0.1 
0.1 
0.1 


3/16 
4/15 
5/18 


6.8 
9.3 
10.76 


17.25/11 
33.1/18.5 
48.5/25 


36.4/25 
80* /55 
90* /80 


68/50 


0.01 
0.01 
0.01 


4/20 
5/18 
6/21 


5.76 
7.3 
8.55 


14.7/9.4 
23.86/14 
33.3/18.5 


32/22 
61/40 
80* /60 


61/45 


5 x 10" 4 


6/28 


5.76 


14.7/9.4 


32/22 


63.5/45 



9^o > 0.03. For fixed temperature correction order, k mlix , one can 
therefore try to find an optimal combination of pivot temperatures 
to cover a large range of temperatures. This is not the absolute min- 
imum with respect to the number of moments, but optimization in 
this way still is beneficial while remaining sufficiently simple. 

In Table Q] we summarize the results of our efforts to cover at 
least the temperature range < kT e < 60 keV for a given precision 
and £ max . We defined different regions of temperatures making sure 
that close to the boundaries the condition on the precision is met 
with some 10% - 20% margin. Far away from the boundaries of 
the different temperature regions the approximations are typically 
much more accurate. The setting for accuracy goal 5 x 10~ 4 A/ S is al- 
ready close to the numerical precision of our approximations and is 
mainly meant to provide an extreme setting for comparisons. Fur- 
thermore, if the electron temperature is smaller than some maximal 
temperature, 7" e max , the total number of required variables can be 
further reduced by dropping moments in regions with T e > r e ,max- 
For a given accuracy goal this defines an optimal value for £ max . We 
also found that the same settings work when < /3 C < 0.01. 

One point we mention is that the settings given in Table Q] are 
in fact independent of the chosen optical depth, t. This means that 
only the scaling of the approximation with electron temperature 
affects the precision. For instance, if the optical depth along a given 
line-of-sight is r - 2 x 10~ 3 , but the temperature is fixed, then 
the absolute precision of the approximation at accuracy goal in the 
second category (denoted with 0.1A/ S ) is actually < 0.02A/ S . 



3 SZ SIGNALS FOR VARIOUS CLUSTER MODELS 

To demonstrate how to use and intetpret the temperature-velocity 
moment method we now discuss the SZ signals for different cluster 
models. We start with the simplest case of an isothermal cluster and 
then work our way through several instructive examples, also intro- 
ducing the simpler moment method that is applicable to sufficiently 
smooth (low temperature-velocity variance) cluster atmospheres. 



3.1 SZ signals for isothermal clusters 

Traditionally, galax y clusters have been modelled usin g a simple 
isothermal /?-model dCavaliere & Fusco-Femianolll978l) : 

N e = /V c .,,[l + (r/r c ) 2 r 3 ' J/2 and T c = const, (14) 

where /V e .o - 10~ 3 cirT 3 is the typical central number density of 
free electrons, r c - 10 kpc is the typica l core radius of clusters, 
and y8 ^ 2/3 (e.g., see iReese et al. In the absence of bulk 

velocities, one therefore has the temperature moments y {k) = #J +1 t 
for r e < r c ,iow and z {k) = N{6 c )(^,t for T e > T siaw . This shows 
that the spatial morphology of the SZ signal is fully determined 
by the overall optical depth factor, r(y) (we sh all ignore small cor- 
rections caused by multiple scattering, e.g. , see lDolgov et alj200ll : 
lltoh et af]|200ll ; IColafrancesco et alj|2003l) . with the same spectral 
shape for any line-of-sight through the cluster. This also implies 
that the average SZ signal measured for an unresolved cluster in 
this case is determined by only one spectral function, and the spa- 
tially averaged optical depth, (t). 

Allowing the cluster to move with a peculiar velocity, /? u 
relative to the CMB one readily obtains the velocity moment^j 



bf = bf 



PcHc^T and bf « /3lP 2 (jJ c )6 k T for clusters 



with T e < r eJow , and cf = fiNme^T, c t 



(*) 



~ PcHcM{e e )^T, and 
c\ k) « y8^P2(/Jc)^V(fe) #e T f° r not clusters. Again the spatial depen- 
dence of the SZ signal factors out and is determined by the one of 
the line-of-sight optical depth alone. It is however clear that the 
SZ morphology becomes frequency-dependent once the cluster 
no longer is isothermal or significant internal motions of the ICM 
are present. The reason for this frequency-dependence is related to 
the variation of the temperature-velocity moments along different 
lines-of-sight, as we illustrate below (Sect. [34). 



3.2 SZ signals for smooth density and temperature profiles 

In more realistic cluster models, the variation of the tempera- 
ture also has to be included. One commo n possibility assumes a 
polytropic temperature profile, 7* e oc p ga 7 dMarkevitch et alj| 19991 : 
iFinoguenov et alj|200ll ; IPratt & Amauj2002l) . Alternatively, one 
can consider fits to the observed t emperature and densi ty profiles 
derived from Chandra X-ray data dVikhlinin et al.l2006|) . 

For our discussion it is only important that the associated pro- 
files are very smooth. This suggests that a good approximation for 
the SZ signal can be found by computing average values for the 
temperature and velocity along the line-of-sight. Corrections to this 
lowest order approximation can then be included using a Taylor- 
series around the average values. Assuming that /? G = 0, we can 
introduce the SZ-weighted electron temperature, 



Te.sztr) = (m e c 2 /k)y m /y { - 1} 



'I 



r e dr. 



(15) 



The integrals are defined by Eq. ([3} but here we do not impose 
any condition on the electron temperature. One can furthermore 
introduce the isothermal temperature moments 



[^r e ,Sz(f)/'«eC 2 ]* +1 T(f). 



(16) 



In general these moments are not identical to y {k) (y) for k > 0, and 



7 We neglec ted tiny corrections caused by deviations from the flat-sky ap- 
proximation dChluba et al.l2012T) . 

8 Here we mean that the spectral dependence of the SZ signal no longer is 
independent of the line-of-sight. 



© 0000 RAS, MNRAS 000, 000-000 



SZ signal processing 1 




x = hv I kT 

Figure 3. Spectral functions Sf? in units of r and So = (2/!/c 2 )(Ar7o//i) 3 » 
270 MJy sr . For comparison we also show the main SZ signal (solid black 
line). The arrow indicates the direction of increasing temperature, with the 
lines within the groups being separated by AT e = lOkeV. 

the ratio p {k \y) = y {k) ly { ^ = T k y {k) /[y {0, ] k+ ' can be used to quan- 
tify departures from isothermality: in regions with p (k> t 1 one 
expects the SZ signal to be poorly represented by just using the 
SZ-weighted electron temperature and line-of-sight optical depth. 
The ratios p fk) are also independent of the overall temperature and 
density scales. They only depend on the shapes of the cluster tem- 
perature and electron density profiles. 

We can write this more formally by using the isothermal mo- 
ment vector, ffij so (r, r e> sz), and expanding the average SZ signal 
around T c sz and t: 

where the first order derivative term canceled after performing the 
line-of-sight average. Here S 1so (t,T cSZ ) = F /n lso (r, 7/ e ,sz) is the 
leading order, average SZ signal, while b* g S K0 (T C ) = Fdjj /Hj so (r c ) 
is the derivatives of the SZ signal with respect to 6 e . To simplify the 
notation, we furthermore introduce 

a/' = U - w +i dr = >{« y (* + 1 ) (-i/ + '-v'"- l) . 

J m=o' m ' 

This means Ay m = 0, Ay m = y (I) - = yP Ap (1) , Ay <2) = 
y^Ap (2) - 3y™Ap (1) , and so on, with Ap w = p {k > - 1. Defining 
_ A/'Vy* = T*A3;*'/ty (0) ]* +1 we then can finally write 

s « s™(t, r c , sz ) + s®(t, r e , sz ) ^ (1 > + s2(t, r e , sz ) ^ (2 > + ... (17) 

where S*(t, T) = (T k /k\)d^Si SO (T, T). In this parametrization the 
observables for the SZ measurement are t, r e _sz, and the tempera- 
ture moments a) {k \ 

With SZpack it is straightforward to compute the required vec- 
tors, SP, with very high precision. In Fig.[3]we show the first few 
S[ '. The typical amplitude of the S[ ' is dropping with k, indi- 
cating that unless the temperature moments a> {k) increase strongly 
with k, higher order terms remain small with the largest contribu- 
tions at high frequencies. This shows that the considered expansion 
becomes fully perturbative unless rather large deviations from the 
smooth temperature profile case are present. 

For example, using the simple fits for one of the hottest clus- 



ters ( T c ^ 9.2 keV at r - 33 kpc), A2029, from the cluster sam- 
ple of IVikhlinin et alj <2006l) . we find D m * 0.16, w <2) a* 0, and 
w (3) ~ o.05 close to the cluster center. This indicates that higher or- 
der corrections decay rapidly. In fact, the correction related to 
contributes at the level of a few percent to the average SZ signal, 
while higher order moments are negligible. We find that even for 
more realistic cases from simulated clusters only a few moments 
of the temperature field need to be known to accurately describe 
the SZ signal (Sect. [3~l4l . Also, in the more extreme case of a two- 
temperature plasma only the first few terms are required (Sect. [3~3l >. 

Some interesting frequencies for the functions S\® are related 
to their nulls, maxima and minima. For SP (which is not shown 
in Fig. [3j we find a rather temperature-independent minimum at 
x ^ 2.26. At x — 4 it crosses zero for T c =s lOkeV, while for 
T c = 50keV the null is located at x m 4.89. Its maximum is located 
at x =* 7.2 for T e = lOkeV and at x =* 10 for T e = 50keV. For 
Sj? (see Fig.O on the other hand we find a maximum at x - 2.12 
and the first null at x s 3.5. The position of the minimum varies 
from x = 5.8 for T, = lOkeV to x =s 7.2 for T s = 50keV. These 
properties might be useful when deciding about the locations of 
frequency channels in future SZ experiments. 

We emphasize that for the computation of S; (t, r c SZ ) and 
pS?'(t, r c ,sz) a large number of temperature terms has to be in- 
cluded. Although in the example given above we only find to 
be significant as additional parameter, this is not equivalent to drop- 
ping higher order temperature terms. This point is very important 
when interpreting future SZ data, since otherwise biased results for 
t, r e ,sz, and u> m are obtained. Similarly, one has to include yS c /j c + 
for the analysis, as we discuss in more detail below. 

We also note that Eq. dl7t is applicable even if the real temper- 
ature distribution is not a smooth function. It is only important that 
the variance of the temperature and thermal pressure remains suffi- 
ciently low to warrant decreasing values of the moments, (J- k \ with 
larger k. This condition is usually fulfilled even in more realistic 
cluster models (see Sect. 13.4^ 

3.2.1 The effect of line-of-sight temperature variance 

Above we showed that the dominant correction to the SZ signal is 
determined by the temperature moment This parameter can 
be interpreted as line-of-sight variance or dispersion of the electron 
temperature but weighted by the optical depth of the scattering vol- 
ume element. One can now address the question of how important 
this term is for the interpretation of the SZ signal. In particular, by 
how much are the deduced best-fit values for r c _sz and r affected if 
the contribution from oi l) is neglected for high-sensitivity, multi- 
frequency SZ measurements. 

We can start by writing the SZ signal for T* sz t r e ,sz as an 
expansion of S(t*, T* sz ) around T e _sz and v. 

S(t\ r* sz ) * S<°>(t, T c , sz ) + S2(t, 7V S z) r + s£>(t, TVsz) ©, 

where = (T* sz - T tSZ )/T tSZ , 0r = (t* - r)/r, and we ne- 
glected higher order terms. To determine the best-fit values for t* 
and r* sz one has to compare to S(t, T e SZ , ~ S[ 0) (t, r e ,sz) + 
S?'(t, r c ,sz) and then minimize the squared difference. The co- 
efficients relating T and to then are only functions of tem- 
perature, and it is straightforward to compute the degeneracy coef- 
ficients a T = -Q T /w (l) and a T = 0/w (1) (Fig.|4j, Both r and r c , sz 
are correlated with aj {[) to a similar degree. The degeneracy is very 
close to unity at low temperatures and only drops to about 1 /2 at 
very high temperatures. We find that or T (T e ) a [l+2.7xl0~ 2 r" 86 ] -1 
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T [keV] x = hvl kT 



Figure 4. Degeneracies of t and T e sz with op) . The expected (biased) best- 
fit values are r* a r[l-a T (r c ) w' 1 '] and T* sz « T c ,sz{l+@T(T e ) w"'] when 
analyzing the SZ signal. 



Figure 5. Second frequency derivative of thSZ distortion in units of r and 
So = (2/i/c 2 )(/:7o//i) 3 m 270 MJy sr~'. The arrow indicates the direction of 
increasing temperature, with the lines being separated by AT C = lOkeV. 



and a T ~ exp (-2.6 x lfr 2 T° M ) match the full numerical result 
with 10% precision. With these expressions we can directly esti- 
mate the expected value for T* sz obtained by computing the best- 
fits to the full SZ signal. 

To determine the degeneracy coefficients we used a very dense 
grid of frequency points in the range x = 0.1 to 30. More real- 
istically far fewer independent frequency bins are available plus 
the signal is averaged over some bandwidth and the beam. Fur- 
thermore, foregrounds and the experimental sensitivity at each fre- 
quency are important. All these aspects affect the degeneracy be- 
tween the SZ parameters, as we explain in more detail below 
(Sect. l5,lT l. Nevertheless, the estimate obtained above gives a rough 
scaling for the importance of line-of-sight temperature variations 
for the interpretation of the SZ measurement. 

3.2.2 The effect of velocity terms on the SZ signal 

Thus far we neglected the effect of bulk velocity on the SZ sig- 
nal. However, the effect of (internal) motions can again be in- 
cluded by expanding the SZ signal around the mean line-of-sight 
values. We first define the two velocity components /3 c j| = /3 cy u c 
and p c± = /3 C Jl - of the moving volume element. Since we 
do not consider the effect of polarization these are the only two 
variables that matter for the SZ signal. At lowest order in /3 C one 
is only sensitive to the line-of-sight SZ weighted averaged veloc- 
ity /?c,u,sz = t -1 J/3 c ,n dr. The associated kSZ signal has a spectral 
shape that is very different from the thSZ. It is therefore rather easy 
to distinguish the two. However, the results for the best-fit values 
of r and r c> sz can be strongly biased even for a small (j8 cj | =* 10~ 3 ) 
line-of-sight velocity if neglected in the analysis (see Sect. 15. lb . 

Looking at the total SZ signal, Eq. (TJ, with 1 = yS^ + f3 2 ci _ 
andyS'/M^c) = P 2 C y ~P 2 ±, it is clear that the SZ signal is only sensi- 
tive to/3 2 ± sz = t' 1 f /3 2 ± dr, but not /3 c ,±.sz directly. For simplicity 
(and because the associated terms are extremely small) we shall 
neglect the variation of temperature-velocity cross terms of order 
O(0l(%,) in the discussion here. Expanding around t, r c sz, Ajl.sz, 
and /J Ci± ,sz = 0, the average SZ signal takes the form 

S * S i0) + S (2 V> + C (1) 6r"> + D a> ^ + E i2) f3 2 c±sz + (18) 

ISO ISO ISO ISO iso^c,±,SZ ' V ' 

where we suppressed the function arguments t, r c .sz, and /J c .||.sz- 



We furthermore defined the signal vectors and dispersion variables, 

c® = (k\m\r' < sz < sz a Ajl , sz s iso (r,r c , sz ,A, ll , sz ,/j c , ± . sz ) 

D» = (fc!)- 1 3* c||sz S lso (r,r c , sz , j e,||.sz,yec,x,sz) 
= (fc!)- 1 a* c±sz S lso (r,r c , sz , J 8 c j t , sz , J 8 c , ± , sz ) 

« w =t- 1 f 08c,i|-yS CJI , sz )* +1 dT 

= (^Vr 1 j (T c - T e ,sz) k V3 c ,\\ - Ajisz) dr, (19) 

where the derivatives are evaluated at y8 c . ± _sz = 0- Since in general 
j8 c can vanish, we did not rescale the velocity derivative terms, and 
correspondingly defined the velocity moments without weighting 
by the mean values. 

As we argue below, for future multi-frequency SZ analysis 
the average SZ signal is described very well by Eq. (1181 . This 
implies that the SZ parameter estimation problem for a single 
lines-of-sight boils down to constraining the set of SZ parame- 
ters p = (r,r c> sz,j6cj|.sz,w (1> ,(T (1> ,/<- (1) ,y6 2 ±sz ). Here we ordered the 
parameters according to their expected importance for the aver- 
age SZ signal. Higher order temperature and velocity moments are 
straightforward to add to the problem and maps of the deduced 
quantities then allow constraining the structure of the cluster. Maps 
of the different components of p therefore provide a compressed 
version of future multi-frequency, high resolution SZ data cubes. 
These new variables then directly allow constraining properties of 
the 3-dimensional cluster's atmosphere. 

3.2.3 Effect of angular resolution on the SZ signal 

When deducing the best-fit parameters for t, r c ,sz> and j8 c ,u,sz it 
is also important to include the effect of angular resolution on the 
measured SZ signal. Gradients of the temperature and velocity field 
within the beam affect these values, as we illustrate here. 

One can model the beam as a simple 2-dimensional Gaussian 
with radius, r res , or a beam variance crjj = ^ 2 cs /[2 ln(2)] « 0.721 r 2 ss . 
With the moment method described above we simply have to re- 
place m with (in) in Eq. l !12t . where (...) denotes the average 
over the beam. For a given electron density, temperature and ve- 
locity field, the associated 3-dimensional integrals are in principle 
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straightforward to compute. However, in detail this can already be 
rather demanding, since several integrals have to be evaluated, re- 
specting the conditions for the different temperature regions. 

On the other hand, assuming rather smooth cluster profiles 
and small beam radius the calculation again simplifies signifi- 
cantly. Let the beam be centered on r = (r v0 , r y o) with central 
values t , r c SZ , and /3 c ,n,sz,o- The average SZ signal around r 
can be computed from Eq. d 1 8t and d!9t by replacing r — > (t), 
Te,sz -» ( T e,sz), and /? c j|, sz -» (/3 c ,ii,sz>- Furthermore, one has to 
substitute f% ±sz — > (/3 2 ± sz) ~ /^e±szo an( ^ com pute the beam av- 
erages of toi (I) , <x (1) , and k ([ \ which read: 



At 



(20) 



« (r>- ] <7Yszr 2 if (T. - (r e , sz )) 2 

(* <1) )*<T>- 1 (J' (fa-Vciw)f<h 

(</») = (r>- ] (Je.sz)- 1 If (T e - (r e , sz » (ft,,, - vSc,d,sz» dr\ . 

These dispersions have two main contributions, one from the varia- 
tion along the line-of-sight and the other from the variation among 
different lines-of-sight inside the beam, which can be separated. For 
instance, by writing T e - <r e-S z> = T t - T c ,sz + TVsz - (r c .sz>, to 
second perturbation order we find 

( w <")^^> + (r e , S z)- 2 ((r; sz )-<r e , sz > 2 ) 

(*«) * + C8 cJ ,, SZ >- 2 ((/3 2 , sz ) " <ft,||,sz> 2 ) (21) 
( £ r (1 ») * ^ + <t-, sz )-' ((r e ,szA,ll> - <re,sz> <ftj,» , 



the filter, to lowest order we have S ~ 5i S0 (Jc) + \d\ S ls0 (x) o~ 2 . In 
Fig.[5]we show d\ S i so (x) for different values of the electron temper- 
ature (/? c = 0). It is clear that the average of the SZ signal over the 
filter also leads to another correction term that affects the deduced 
values for t, 7" c ,sz, and /} c j|. The significance of the introduced bi- 
ases depends on the number and position of channels as well as the 
band-width. At low frequencies the effect is very small, but close 
to x ~ 7 it can reach ai 10% of the main signal for band- width 



• SkeV and/ =0.2 



•- A= 1 (ca < "-0.1l,ra <2 '- 0.056) 

— A=2 (0)'" - 0.33, a' 2 ' - 0.28) 

— A = 3 (CO™ - 0.56, a' 2 ' ~ 0.63) 

— A = 4 (0) ll, -0.79, (i)' 2 '- 1.05) 

1 



Figure 6. Importance of higher order moments for a two-temperature 
plasma. We assumed T c i = 5keV and f T = 0.2. The absolute difference of 
the full two-temperature distortion, A/ lwo _7, and the approximation A/ ;lppl 
which only uses the average SZ temperature and the temperature dispersion, 



df l) , is shown in units of A/ s ^0.013 MJy sr 1 (i.e., AX = A/ two _3 
The temperature difference was varied in each case as annotated. 



-A4ppr). 



where, for example, af Q = ay- '(r ), is the line-of-sight tempera- 
ture variance at the beam center. Performing a Taylor expansion of 
quantities around ro up to second order we have 

(x)*x + ^o-l[dl + dl]x 

(X 2 ) - (X) 1 « a\ [(d rx X f + (d,,X ) 2 ] (22) 
(XY) - (X) (¥) * a\ [(d ri X Xd r Jo) + (d rv X )(d ry Y )] 

As these expressions show, the beam average values are affected 
only by terms related to the second derivatives, while the variances 
depend on first derivative terms. Depending on the position of the 
beam with respect to the cluster center this introduces an interesting 
spatial dependence of , (o""'}, and ("■" , ). The precise values 
depend on the structure of the cluster, and only in the limit of very 
small beam radius will they be dominated by the line-of-sight vari- 
ances. When deducing cluster parameters from SZ measurements 
this aspect has to be kept in mind. If this effect is ignored it can 
again lead to biases in the inferred parameters. 

3.2.4 Effect of frequency resolution on the SZ signal 

To estimate the effect of frequency resolution on the SZ signal we 
can use the fact that the SZ signal is a very smooth function of x. 

Defining the average frequency x and dispersion a\ = x 1 — x 2 over one finds oj {[) < u) {2) 



comparable to 10% (assuming a simple top-hat filter). As Fig. [5] 
shows, the other two regions with large curvature in the SZ signal 
are around x = 3 and x — 10-15, with the latter having a significant 
dependence on electron temperature. 

3.3 SZ signal for a two-temperature plasma 

In the previous section we illustrated how the SZ signal is affected 
by the temperature profile of the cluster. However, for the analy- 
sis we used very smooth average profiles which allowed us to ob- 
tain an accurate description of the SZ signal with only a few new 
parameters (see Sect. 13. 2t . In general the variance of the tempera- 
ture, density and velocity field along different lines-of-sight could 
be much larger than suggested by the simple average profiles. 

One illustrative example is to consider a two-temperature 
plasma, with one low temperature, T e> \ = T c , and a high temper- 
ature region at T e 2 = T e (l + A), where A is the relative difference 
between the two temperatures. We furthermore have to specify the 
optical depths of each region, T\ = t(1 - f T ) and r 2 = f T r, where r 
is the total optical depth along the line-of-sight. From this it follows 



r ejSZ = r 6 (i+/ T A) 

jT (l+f T A) 2 Jt 

(l+/ r A) 3 yjT 



(23) 



!)(/;' - 2) 



The temperature dispersion has a maximum at f T = [2 + A] 1 with 
«£L = jA 2 /[l + A]. Furthermore, for A > [1 - and f T < 1/3 



In Fig.[6]we illustrate the importance of higher order tempera- 
ture moments for a two-temperature plasma with T sA and f T = 0.2 
and different values of A. For rather small temperature ratio the av- 
erage SZ signal is approximated very well using the average SZ 
temperature and with higher order temperature moments lead- 
ing to very small corrections. For the most extreme case shown we 
have T c | = 5keV and T] = 0.8 with a high temperature compo- 
nent r e .i = 25 keV and T2 = 0.2. In this case higher order moments 
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Figure 7. First few moments of the temperature field for a simulated cluster from lNagai et al.l l l2007blal) . We only show one projection; r e ,sz is in keV. We 
marked several regions for which we also computed the average SZ signal (see discussion in main text). The highest density peak (marked with an open circle) 
defines the cluster center. The box is 1 Mpc Ir 1 on the side and the cluster mass is as 3.8 X 10 14 M Q . 



Table 2. Optical depth, SZ temperature, and first two line-of-sight temper- 
ature moments for the regions marked in Fig. [7] 



Region 


r[l(r 3 ] 


Te.sz [keV] 




o/ 2 > 





7.7563 


4.1051 


0.066421 


0.021254 


I 


3.1682 


6.3163 


0.14538 


-0.021378 


II 


0.62251 


1.1251 


0.69288 


1.1797 


III 


1.3738 


3.8462 


0.29556 


0.17844 



are important at the level of a few percent for the SZ signal, but 
we found that the approximation improved dramatically when only 



a case that 



adding u> {1) . In this example one also had a> (l 
was not found for the smooth cluster profiles discussed in Sect. 13.21 
Encountering <±> a> < u) {2) therefore can be valued as indication of 
non-smooth, high variance plasma, possibly with more than one 
temperature component. However, to assess the statistical expecta- 
tion for such cases detailed simulations of clusters within the cos- 
mological context should be performed. 



3.4 Thermal SZ signal for a simulated cluster 

We close our discussion of approximations of the SZ signal with 
one exam ple taken from hydro dynamic simulations presented pre- 
viously in lNagai et al.l d2007bl l5). The simulation is performed with 
the Adaptive Refinement Tree (ART) jV -body+gas-dynamics code 
dKravtsov et alj2002l : lRudd et alll2008h . which is an Eulerian code 
that uses adaptive refinement in space and time, and non-adaptive 
refinement in mass to achieve the dynamic range necessary to re- 
solve the cores of halos formed in self-consistent cosmological 
simulations. The simulation assumes a flat ACDM model: Q m = 
1 - O a = 0.3, a b = 0.04286, h = 0.7 and cr 8 = 0.9, where the Hub- 
ble constant is defined as 100/i km s -1 Mpc -1 , and <r 8 is the mass 
variance within spheres of radius 8 h' 1 Mpc. The simulation shown 
here is of a cluster with mass M500 ^ 3.8 x 10 14 M G undergoing 
a near 1 : 1 merger at z = 0.25 jNelson et al ] |2012l) . The simulation 
was run using non-radiative gas dynamics on a uniform 128 3 grid 
with 8 levels of mesh refinement, corresponding to peak spatial res- 
olution of 7 /T 1 kpc. The dark matter (DM) particle mass in the 
region around the cluster was m p 9.1 x 10 8 h' 1 M Q . 

In Figure|7]we show projections of t, r cS z, w (1> , and a> i2) . All 
variables are rather smooth close to the cluster center, but show 
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Figure 8. Dependence of the crossover frequency on 7" c sz and /? c .||.sz- 

significant variations in the outskirts. We marked regions with tem- 
perature variance ^ 0.1 - 0.2 (see Table [2] for details). In 
these regions also cl> (2> acquired significant values, so that one could 
expect large contributions to the average SZ signal from higher 
order terms. The highest temperature region (marked with T) is 
slightly displaced from the maximum of the line-of-sight optical 
depth (marked by an open circle). It is also coinciding with a region 
of enhanced temperature variance. Region 'IF on the other hand is 
one of the coolest regions in the cluster but very large line-of-sight 
temperature variation. 

We computed the SZ signal for the different regions using di- 
rect integration of the distortion for each frequency. We also con- 
firmed that with the moment method developed in Sect.|2]gives the 
same result but at a significantly lower computational cost. In both 
cases, bulk velocity terms were omitted. We then compared the ap- 
proximate SZ signal obtained with the moments summarized in Ta- 
bleland found that the deviations were very small, consistent with 
higher order temperature moments being negligible. This indicates 
that even for more realistic cluster models the SZ signal is well ap- 
proximated using the smooth profile expansion, Eq. ( 118) . However, 
the considered cluster generally has rather low temperature so that 
relativistic corrections are not expected to be as important. We will 
investigate more massive clusters in our future work. 



4 EFFECT OF DIFFERENT LINE-OF-SIGHT 
VARIATIONS ON THE SZ NULL 

One of the important bands for SZ observations is defined by the SZ 
null. It is well-known that the position of the so called crossover fre- 
quency depends on the temperature of the cluster gas and average 
line-of-sight velocity. In particular, higher order temperature cor - 
rections affect its precise location (e.g., see lNozawaetaLlll998bl) . 
This is illustrated in Fig. [8] where the position of the null was com- 
puted using one of the functions implemented in SZpack. The re- 
su lts are in good agreem ent with the simple fitting formula given 
bv lNozawa etai] Jl998bh . 

However, it is clear that the line-of-sight temperature and ve- 
locity variations also affect the position of the null (see Fig. |9). 
Generally, the change in the position of the SZ null is rather small 
in the cases shown. The shift is nearly temperature-independent for 
cr (1) , while for it grows with T c . We also present the effect of 




r [keV] 



Figure 9. Dependence of the crossover frequency on T^sz, o>^\ cr' 1 ', ic \ 
and j8^ ± sz - We present the absolute shift, x c - v c o, where x c o is the case 
with/3 cJ | = 6J (1) = tr (1) = ff (1) =/3j? ±sz = 0. 

second order velocity terms on the crossover frequency. The line- 
of-sight variance of the velocity causes a decrease of the crossover 
frequency, while yS^ ± sz leads to an increase. Both effects are rather 
small for a typical cluster temperature, even for the relatively large 
fiducial values that were chosen in the illustration. They further- 
more practically cancel each other, indicating that for SZ parame- 
ter estimation they are degenerate. In the simulated cluster shown 
in Fig.|7]we found maximal values of 01 sz ^ 10~ 5 and/e (1) 10~ 5 . 
However, this depends strongly on the chosen projection and post- 
merger time. It is therefore important to analyze a larger number of 
simulated clusters to estimate the typical values for 01 L sz and 

4.1 Effect on the deduced SZ temperature 

The results of Fig. [9] indicate that conclusions on the temperature 
of the cluster drawn from measurements close to the SZ null can be 
affected by line-of-sight temperature and velocity variations. For 
instance, assuming /? c n = 0, we find 

x c ~ 3.83[1 + 0.01 lf c -8.6 x 10^ 5 f c 2 ] 

Ax c a 0.042f c [l - 0.069f c ] (24) 

with f c = r c /5keV for T e < 30keV. This means that for SZ 
temperature T c .sz = 5 keV and =2 0.2, we find x c ~ 3.88, 
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which would imply a temperature of r eSZ m 6keV when assum- 
ing w (1> = in the analysis. This is a =i 19% bias towards higher 
temperature. More generally we find a bias of 



1 - 0.0697; 
1 - 0.016f,, 



„(i) 



(25) 



for T e < 30keV. The degeneracy between T e and <x) (1) is very close 
to unity at low temperatures; however, for r e .sz = 30keV and 

= 0.2 the temperature bias reduces to ^ 13%, which is also 
in good agreement with the behaviour found in Fig. [4] Although 
the absolute shift in jc c caused by a> {1) increases with temperature 
(cf. Fig. [9}, the relative shift decreases. 

One can also estimate the effect of higher order velocity terms 
on the inferred value of r c _sz- For contributions related to <9(/3 2 ) the 
effect is very small and only relevant for very low temperature gas, 
if at all. Similarly, the value of a m is expected to be much smaller 
than < 10~ 3 , unless high temperature gas also has large bulk ve- 
locity. For example, for the simulated cluster shown in Fig. [7] we 
found maximal values of tr (1 ' ^ 4.7 x 10~ 4 . Therefore the correc- 
tion caused by is expected to be dominant, unless the overall 
temperature of the electron gas is very low. 

As explained in Sect. 13.231 also variations caused by averag- 
ing the SZ signal over the beam introduce a dispersion. This has the 
same effect as the line-of-sight temperature variation, but for small 
beam the latter should again dominate. 



4.2 Effect of frequency resolution on the SZ null 

It is also straightforward to estimate the effect of the bandpass on 
the SZ null. The shift in its position is approximately given by 



Ax. 



~2d A ,S(jt c ) * 

-0.2[1 - 0.32f c + 0.03lf c 2 - 0.0018f e 3 ] <t\ 



(26) 



for T e < 30keV. At T e = 23keV the shift in the crossover fre- 
quency caused by the average over the filter becomes small. Also, 
for a top-hat frequency filter one has <x 2 . = Ax 2 /12 1.22(Av/v c ) 2 , 
which for a bandwidth of 10% and T c ^ 5 keV implies a shift of 
Ax c ^ -0.002 in the crossover frequency. This is expected to result 
in a a 4% bias of T c towards lower values. More generally we find 



AT e 

17 



5.8 /Av\ 2 1 - 0.298f c + 0.0247; - O.OOllf 3 



1 - 0.0167; 



(27) 



caused by the band-width of the frequency channel around x c (r c ). 

One should mention that we assumed that the frequency band 
is centered exactly on the crossover frequency, which of course 
changes with temperature, T c . If the central frequency differs from 
x c the response is also affected and hence the deduced temperature 
r c . In the limit of very narrow bandpass this problem disappears, 
but otherwise this effect should be taken into account when inter- 
preting SZ data at high precision. 



5 SZ PARAMETER ESTIMATION 

In Sec.[2] we developed a new method for accurately calculating the 
mean SZ signal with relativistic corrections, and in Sec. [3] we ex- 
tended this to include fluctuations in temperature and velocity along 
the line-of-sight for smooth cluster profiles. Subsets of the model 




parameters {t, T e ,sz, oi 



,(1,2,3) ,-(1,2,3) „(1,2,3) 



ft,||,SZ>A. i < 



can then 
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Figure 10. Joint marginalized posterior distribution for a two-temperature 
model of the Bullet cluster showing the fractional optical depth f T which 
can occupy a component of gas with temperature AT higher than the main 
component. Here, the temperature of the main component is (7.2±3.9) keV. 
When the additional hot component is 20 keV hotter, it can occupy a signif- 
icant fraction of the optical depth while for larger excursions in temperature 
is must have progressively smaller optical depth. Note that non-zero AT is 
driven by the measurement at 857 GHz, where point source contamination 
could be significant. 

be estimated from measurements of the SZ spectrum^ Here we 
describe an estimation of these parameters and their interpretation 
using a compilation of data and simulatio ns. The parameters are 
estim ated in a Markov Chain Monte Carlo dForeman-Mackev et alj 
2012) which is available with SZpack and can be adapted for par- 
ticular experimental settings. 

Rather than performing an exhaustive survey of experimental 
configurations and considerations, we describe a few informative 
cases. Even for a given combination of multi-band data, a com- 
plete interpretation is beyond the scope of this paper because it 
demands a consideration of contaminants, calibration errors, and 
offsets which can all be correlated across the bands. Particulars 
of the ob servation strategy m ay also need to be fully modelled 
(e.g., see IZemcov et alj|2012l ). The examples here are therefore 
meant to show some degeneracies and interpretative issues that 
arise in multi-band SZ parameter estimation. The same conclusions 
broadly apply to stacked cluster studies, except there, for example, 
the kSZ may average to zero. 

5.1 Interpreting temperature variance in a recent SZ data 
compilation for the Bullet cluster 

The simplest estimation problem which demonstrates evidence of 
temperature dispersion is the two- temperature model of Sec. 13.31 
Here the parameters are a reference temperature T e , a temperature 
difference AT, and the fractional optical depth f T < 1 in the region 

9 We will drop the subscript 'SZ' for convenience, but it is important to 
bear in mind that for instance the derived electron temperature, T c = T c sz, 
i.e. it is weighted by the electron number density. 
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with temperat ure T e +AT. A two-temperatu re case was also recently 
considered by IColafrancesco et all d20 1 lh as model for the Bullet 
cluster gas. The two SZ contributions from the gas at T c and T C +AT 
are computed to all relevant relativi stic correction orders using 
SZpac k. Following the compilation of IProkhorov & Colafrancescol 
ilOlA we consider data across the bands ( 150, 275, 600, 857) GHz 
with mean values and errors of A/ =(-0.325 ± 0.015,0.21 ± 
0.077, 0.268 ± 0.03 1, 0.0 97 ± 0.019} MJ y /sr fro m observations of 
iGomez et al] J2004) and IZemcov et ail J2010h. We assume that 
t = 0.0138 ± 0.0016 is known ( IProkhorov & Colafrancescoll2012h 
as a Gaussian prior probability distribution from X-ray observa- 
tions. This helps break the strong degeneracy between t and T s . 
For the latter we find T c = (7.2 ± 3.9) keV. Figure [10] shows 
the marginalized posterior distribution for AT = 25+ j' keV and 
f T = 0.27+jj 2 * (where two-sided errors are given for non-Gaussian 
distributions at percentiles (16%, 50%, 84%), equivalent to l<x in a 
normal distribution). Conclusions from a model with yS c ,n (taken to 
be the same for both components) as a free parameter are similar at 
this level of uncertainty. 

We next consider an ICM model with line-of-sight tem- 
perature dispersion about the mean using the free parameters 
(t, r e ,w (1) ), where cm {[) = (T 2 )/[tT 2 sz ] - 1 is a measure of the 
temperature dispersion (see Sect. 13.2) . In this case, we obtain 
T e = (13.5 ± 1.8) keV and co il) = 0.8 ± 0.4 with a modest corre- 
lation corr(r c , = -0.4. When yS c || is added as a free param- 
eter, we find T e = (13.1 ± 2.0) keV and w m = 1.0 ± 0.6 with 



(-1.7 ± 2.1) x 10 3 and corr(/? L ,||, <o w ) 



-0.6, showing 



that there is no significant indication of non-zero line-of-sight bulk 
velocity. We can extend the fit to include or 2 ', a measure of the 
skewness in the temperature distribution. Here, to {>> =1.1 ±0.6 and 
u> {2) = 2.0 ± 2.6 (no strong evidence for more than dispersion) with 
a strong correlation, con^o/ 1 ', a> <2) ) = 0.75. It is also informative to 
compare with a model where <y (1) is fixed to zero and only T e and 
t vary. In this case, T c = (14.8 ± 2.1)keV, but as shown in Fig. II II 
this is a poorer fit to the data because it lacks the freedom to explain 
the flux measurement at 857 GHz. 

The two previous models depart from the first reported mea- 
surem ent of temperature dispersion bv lProkhorov & Colafrancescol 
< |2012|) , where the temperature expansion is used directly to con- 
strain the variance: 



AI 

J- ~ T\g (x) 



(kT e ) {(kT e ) 2 ) <(fcr„) 3 ) 
2 + gi(x) — TT~ + 8 l(x) — TT~ 



(28) 



Here go, gi and g2 are the spectral functions of the asymptotic ex- 
pansion given in Ichallinor & Lasenbvl Jl998h and flux normaliza- 
tion, /„ = (2h/c 2 )(kT /h) 3 x 270MJy sr~'. The temperature vari- 
ance can be inferred directly through cr 2 = ((kT s ) 2 ) - (kT s ) 2 , by fit- 
ting for the coefficients of go, gi and g2. As a check of the MCMC 
here, we find a temperature variance of (9.2 ± 2.5) keV, in agree- 
ment w ith (9.5 ± 2.6) keV reported in IProkhorov & Colafrancescol 
The average temperature of the Bull et cluster inferred from 
X-ray observations dMillion & Alien! l2009l) is 14.5 keV. At these 
temperatures, Fig.[2]shows that the asymptotic expansion is a poor 
approximation even to order {(kT e ) 10 ) I (m e c 2 ) i0 , which worsens for 
higher temperature excursions. Because the first few moments do 
not adequately describe the mean relativistic corrections, it is also 
difficult to interpret the dispersion inferred from the first two mo- 
ments in Eq. J28l >. Nevertheless, converting the moment constraints, 
we find that T, = (13.7 ± 1.9) keV and ~ a- 2 /T 2 = (0.5 ± 0.3) 
for this model, which is consistent with the dispersion measured 
above. However, as the experimental sensitivity and frequency cov- 




Figure 11. Mean SZ spectra for four models of the data (black points) 
I Prokhorov & Colafrancesco 2012). A model including only 7" c with rela- 
tivistic corrections at all relevant orders (yellow) does not have the freedom 
to describe the 857 GHz flux, and drives the temperature up. An asymp- 
totic expansion of relativistic corrections to order ((&r c ) 3 )/(m e c 2 ) 3 (red) 
describes the data, but is insufficiently converged to self-consistently ex- 
plain the relativistic corrections to the SZ flux. The green curve is a model 
with two regions of different temperature and the blue curve is a model 
with dispersion about a mean temperature (both cases include all relevant 
orders in the relativistic corrections). Measurements around and above the 
SZ maximum have the potential to differentiate these models, while close 
to the SZ null and below all models practically coincide. 



erage increases the model given by Eq. Q28J should lead to biased 
values for T e and a/ l \ 

In any of the models above, the detection of dispersion is 
driven by the flux measurement at 857 GHz. Removing that mea- 
surement from the two-temperature model parameter estimation re- 
sults in an upper bound rather than detection (AT < 56 keV, 95% 
confidence). In the dispersion model with parameters \T c ,t, w (1> ), 
we find w ( " = 0.4 ± 0.7. Most of the constraint on relativistic ef- 
fects and dispersion comes from above the SZ null, but in this fre- 
quency region it becomes more difficult to separate point source 
contamination, both because IR fluxes increase at higher frequency 
and because both the SZ and contamination are positive fluxes. 

The parameter estimates here are based on flux measurements 
in four frequencies, and yield qualitatively similar conclusions for 
the current sensitivity. Across a range of frequencies, however, the 
three models do predict different SZ spectra, shown in Fig. Qj] 
Here, the expression of Eq. ((28) shows the most variation with re- 
spect to other models with temperature dispersion because three 
temperature moments in the asymptotic series do not fully describe 
the relativistic corrections and temperature dispersion. The other 
two models show slight differences above the SZ maximum, in- 
dicating that more precise measurements at a wider range of fre- 
quencies could help differentiate between physical interpretations 
of the dispersion. Still the differences in the SZ signal above the 
null are not much larger than ^ 5% - 10%. There are limited trans- 
parency windows in the atmosphere, but at frequencies above the 
null, instruments can achieve similar resolution and sensitivity with 
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comparatively smaller apertures, so that balloon and space-borne 
instruments could have the potential to constrain these frequency 
regions. Note that the distortion around the null is virtually identi- 
cal in all the models described here. 

In the constraints above, we neglected the finite bandwidth of 
the observations. It is customary to define an effective band center 
for a given spectral index as the frequency at which the measured 
flux matches the true flux. We can expand smooth spectra about 
this effective band center and see that the linear term drops out and 
that the change in center is due to the second derivative curvature 
of the spectrum across the band (see Sect 14. 2t . The extent of this 
effect depends on experiment-dependent parameters, but as a rough 
estimate, a 10% bandwidth near the maximum at x ~ 7 can produce 
a nearly 10% bias. In the data considered above, for a 10% band- 
width, we find T c = (13.6 + 1.8) keV and a) (i > = 0.72 + 0.4, a shift 
which can be neglected here, but should become important as the 
measurements improve. 

We can also examine the impact of improved measure- 
ments on the inference of Fig. [12] considers the param- 
eter set (t, T e ,tj (1) ) under improv ements to the dataset used in 
IProkhorov & Colafrancescol d2012l) . By reducing the SZ measure- 
ment errors a factor of five for the same r prior, the a> il) - T e cor- 
relation increases to -0.7 and the parameter errors improve mod- 
estly. Both the correlation and residual uncertainty are driven by 
the imperfect r prior, and deeper SZ measurements have little ad- 
ditional benefit. Shrinking the t prior by a factor of five on top of 
the improved SZ flux measurements significantly improves these 
constraints. This also emphasizes the importance of combining 
SZ measurement with X-ray observations to break degeneracies 
among model parameters. 





■%■ 


/try j 


0.0 0.5 1.0 1.5 












B 10 12 14 16 IB 20 22 

T e keV 





Figure 12. Posterior distributions of {r e , r, fa/ 11 ) achieved by improve ments 
to the measurements compiled in lProkhorov & C olafrancesco (201^). The 
green solid line gives ID marginalized probabilities and 2D contours (95% 
enclosed) for the original dataset. This shows a clear indication of scatter 
about the central temperature through significantly non-zero Improv- 
ing the SZ flux measurements by a factor of five (keeping the measured 
central values) only improves the constraint on T c and (J- > (blue dashed 
line) modestly because of their degeneracy with t. Additionally shrinking 
the t prior errors by a factor of five (red dot-dashed line) breaks this to yield 
the strongest constraints. 



6 DISCUSSION AND CONCLUSION 

We developed a novel method for extracting information about 
the state of the ICM from SZ observations which is based on 
moments of the cluster's temperature and velocity field. Both a 
quasi-exact approach (Sect. |2), which should be useful for com- 
puting the precise SZ signal from simulated clusters in an effi- 
cient manner, and an approximate treatment based on the smooth- 
ness of the cluster profile (Sect. 13. 2t were considered. In fu- 
ture SZ parameter estimation and cluster profile reconstruction 
the latter method should be widely applicable, with the main SZ 
observables consisting of a small number of parameters, p = 
(t, r c _ sz , cD il23 \ cr (I - 2 - 3) , K (1A3) ,^ c ,[],sz,^ ± sz'- Ma P s of tnese param- 
eters allow very precise representation of the SZ signal and hence 
a compression of the observational data to only spatially dependent 
variables. These can be directly used in the inversion problem that 
determines the cluster temperature and velocity structure. 

Our discussion of different cluster models indicated that the 
above compression of the data should work even in more general 
cases, when the line-of-sight temperature and velocity distributions 
become non-trivial. The reason is simply that the SZ signal itself is 
very smooth in frequency and higher order temperature and veloc- 
ity moments contribute in a diminishing manner, once the average 
SZ signal (which includes a large number of temperature correction 
terms) is taken out. We also showed how the position of the SZ null 
is affected by temperature variations along the line-of-sight. If ne- 
glected this can lead to biases in the deduced electron temperature 
at the level of =; 10% -20% (see Eq.|25t. Velocity correction terms, 
beyond the normal kSZ effect, are generally expected to be less im- 
portant; however, a more detailed study using hydrodynamic cluster 



simulations is needed to give a more quantitative answer about the 
typical magnitude of different terms. 

We furthermore illustrated the importance of high frequency 
channels for distinguishing different c luster atmospheres (Fi g. [Tfl 
using recent data for the Bullet cluster jProkhorov & Colafrancescol 
|2012|) . Our analysis indicates that presently at two-temperature 
model remains indistinguishable from a cluster model with line- 
of-sight temperature variance, but a single temperature model is in 
tension with the measured SZ flux at 857 GHz. We also explain 
why the results for the temperature dispersion deduced with an 
asymptotic expansion of the SZ signal are difficult to interpret. In 
the future, it will be important to consider foregrounds and sys- 
tematic effects in more det ail. Also additional physical effect, e.g., 
multiple-scattering effects fc>olgovetai1l200lh or the SZ effect for 
non-thermal electron populations l Enfilin & Kaisedbood) . should 
be taken into account. We plan to extend SZpack with these fea- 
tures in another publication. 

To realize the full power of SZ measurement as a cosmolog- 
ical probe, studies of individual clusters and their atmosphere will 
become very important. It is especially critical to extend the studies 
of thermodynamic and velocity structures to the outskirts of high- 
redshift clusters in order to understand the biases in the global prop- 
erties of galaxy clusters and their evolution across the cosmic time. 
Because the SZ signal is independent of redshift and linearly pro- 
portional to the gas density, high-resolution SZ observations can 
enable detailed studies of the ICM structure in the outskirts of 
high-redshift clusters, and hence are highly complimentary to X- 
ray observations whose sensitivities are limited to the inner regions 
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of clusters or the outskirts of nearby systems. The combination of 
future spectroscopic X-ray measurements with high resolution SZ 
data will open new opportunities for studying cluster astrophysics, 
and the moment method developed here provides one important 
step in this direction. 
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APPENDIX A: DEFINITION OF FUNCTIONS FOR 
ASYMPTOTIC EXPANSION 

From Eq. (25) of CNSN we can directly identify 



M iow _ y*m + z) xYk - -Dim 

6 3 
D [ ° w = Df n -xd x Y k 



(Al) 



The functions Y k , K t kln , flf 1 , and Qf n are all defined in CNSN. The 
differential operators are D x = 3xd x + x 2 5 2 , T)\ = 2 + xd x , and 
D" = 1 - xd x . With the relations given in CNSN it is straightfor- 
ward to show that 
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fc=0 
2n+2 



= ~Yj 4° - d A * + 2k A i +1 + A i +2 ] «k w 

with n Pl (x) = - 1] and A* = x k d k . This implies 
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mi™ = i Yjtfk - 4 n,) i k(k + 2) A * + {2k + 3) A * +1 + A * +2 ] 
D T = YM? - 4"') [ k K + A f '] «hW 



(A3) 



k=0 
2n+2 



QT = t 2 (a ?' + <?f - 2 4 n) ) [*(* - u a* + 2*A* + 1 + a a ; +2 ]« p1 



We extended the sums for Y k and its derivatives to k = using 
af = for n > 0. 



APPENDIX B: DEFINITION OF FUNCTIONS FOR THE 
IMPROVED BASIS SET 

It is straightforward to transform the set of basis functions given in 
CNSN from the cluster frame into the CMB rest frame. However, 
to obtain the expressions for Y^' gh , M^' gh , D^ gh , and 2j' sh , a few 
intermediate steps are needed. Using the definitions of CNSN gives 

Y k ~ £ k 



= cf n -xd x z k 



(Bl) 



Q 



high.. = s m + i 2 d 2 Zk + Izrcf 1 . 



These can be used like the functions Y^' gh , M^' gh , Z)J' gh , and Q h k ' gh in 
Eq. but with N(8 C ) &\ -> N(Q S ) (0 e - 6,$) k I N{8 tfi ). Since 

k 

"V-Befif^NmV! (B2) 



-y( k 



^)Q(-ft.or 



by defining the lower triangular matrix, Tjj 
we have the transformation X h ' gh = T T X h ' gh -* to the basis Y^ gh , 
Mf gh , D^ gh , and Q^ g \ Here we used (X^*f = (< gh '\ ...,X^) 

and (Z hlgh ) r = (< sh , ...,<t), with X 6 {7, Af, O, Q}. 

Finally, with the definitions of CNSN, the integrals and deriva- 
tives in Eq. l lB 1 1 can be written as 

dV dV A ^ 



M h : g \x) : 



N_r 

3k\J 



dju d// d k ff t 



— [At(x') - M(x)] dn&n'rf-da], 



k k\J d/jdfi' 



2_ gke-Ay/ei 



gW&ii&ti'rf&ri 



(B3) 



N C dVp 

_/v r dV 2 

" 3Wj " 

yv r 



[^(x')-^(x)]d^d/77 2 d77, 



dju d//' d k ff t 
d 2 o-o 

d// dyu' d k ff- e 
2N C d 2 cr! 3 ,: e- A '' /fl « 



"f- 

t!j c 



— [Q(x') - A Q(x)] d^d^'^ 2 d^ 
[(3(x')-QW]d^d^'^ 2 d77 
CJ(x')dyudyu'^ 2 d^. 



3fc!j d^d/i' 3*6^ 
We implemented these alternative basis functions for SZpack. 

REFERENCES 

Benson B. A.. Church S. E., Ade P. A. R., Bock J. J., Ganga K. M., Henson 

C. N., Thompson K. L., 2004, ApJ, 617, 829 
Benson B. A. et al., 2003, ApJ, 592, 674 
Birkinshaw M., 1999, Physics Reports, 310, 97 
Birkinshaw M, Gull S. E, Hardebeck H., 1984, Nature, 309, 34 
Birkinshaw M., Hughes J. P., Arnaud K. A., 1991, ApJ, 379, 466 
Carlstrom J. E., Holder G. P., Reese E. D., 2002, ARA&A, 40, 643 
Cavaliere A., Fusco-Femiano R., 1978, A&A, 70, 677 
Challinor A., Lasenby A., 1998, ApJ, 499, 1 
Challinor A., Lasenby A., 1999, ApJ, 510, 930 
Chluba J., 201 1, MNRAS, 415, 3227 
Chluba J., Hiitsi G, Sunyaev R. A., 2005, A&A, 434, 81 1 
Chluba J., Mannheim K., 2002, A&A, 396, 419 
Chluba J., Nagai D., Sazonov S., Nelson K., 2012, MNRAS, 426, 510 



© 0000 RAS, MNRAS 000, 000-000 



16 Chluba et al. 

Colafrancesco S., Marchegiani P., Buonanno R., 2011, A&A, 527, LI 

Colafrancesco S., Marchegiani P., Palladino E., 2003, A&A, 397, 27 

Diego J. M., Mazzotta P., Silk J., 2003, ApJL, 597, LI 

Dolgov A. D., Hansen S. H., Pastor S., Semikoz D. V., 2001, ApJ, 554, 74 

EnBlin T. A., Kaiser C. R., 2000, A&A, 360, 417 

Finoguenov A., Reiprich T. H., Bohringer H., 2001, A&A, 368, 749 

Fixsen D. J., Cheng E. S., Gales J. M., Mather J. C, Shafer R. A., Wright 

E. L., 1996, ApJ, 473, 576 
Fixsen D. J., Mather J. C, 2002, ApJ, 581, 817 

Foreman-Mackey D., Hogg D. W., Lang D., Goodman J., 2012, 
ArXiv: 1202.3665 

Gomez P. et al., 2004, in American Institute of Physics Conference Series, 

Vol. 703, Plasmas in the Laboratory and in the Universe: New Insights 

and New Challenges, Benin G., Farina D., Pozzoli R., eds., pp. 361-366 
Hand N. et al., 2012, Physical Review Letters, 109, 041 101 
Hansen S. H., Pastor S., Semikoz D. V., 2002, ApJL, 573, L69 
Hughes J. P., Birkinshaw M., 1998, ApJ, 501, 1 

Itoh N, Kawana Y., Nozawa S., Kohyama Y, 2001, MNRAS, 327, 567 
Itoh N, Kohyama Y, Nozawa S., 1998, ApJ, 502, 7 
Itoh N, Nozawa S., 2004, A&A, 417, 827 

Kitayama T., Komatsu E., Ota N., Kuwabara T., Suto Y, Yoshikawa K, 

Hattori M., Matsuo H., 2004, PASJ, 56, 17 
Komatsu E. et al., 2001, PASJ, 53, 57 
Komatsu E., Seljak U., 2002, MNRAS, 336, 1256 
Korngut P. M. et al., 201 1, ApJ, 734, 10 
Kravtsov A. V., Klypin A., Hoffman Y, 2002, ApJ, 571, 563 
Lamarre J. M. et al., 1998, ApJL, 507, L5 
Markevitch M., Vikhlinin A., 2007, Physics Reports, 443, 1 
Markevitch M., Vikhlinin A., Forman W. R., Sarazin C. L., 1999, ApJ, 

527, 545 

Marriage T. A. et al., 2011, ApJ, 737, 61 

Menanteau F. et al., 2012, ApJ, 748, 7 

Million E. T., Allen S. W., 2009, MNRAS, 399, 1307 

Mroczkowski T. et al., 2012, ApJ, 761, 47 

Nagai D., Kravtsov A. V, Kosowsky A., 2003, ApJ, 587, 524 

Nagai D., Kravtsov A. V, Vikhlinin A., 2007a, ApJ, 668, 1 

Nagai D., Vikhlinin A., Kravtsov A. V, 2007b, ApJ, 655, 98 

Nelson K., Rudd D. H., Shaw L., Nagai D., 2012, ApJ, 751, 121 

Nozawa S., Itoh N., Kawana Y, Kohyama Y, 2000, ApJ, 536, 31 

Nozawa S., Itoh N., Kohyama Y, 1998a, ApJ, 508, 17 

Nozawa S., Itoh N., Kohyama Y, 1998b, ApJ, 507, 530 

Nozawa S., Kohyama Y, Itoh N, 2009, Phys.Rev.D, 79, 123007 

Planck Collaboration et al., 201 1, A&A, 536, A8 

Pointecouteau E., Giard M., Barret D., 1998, A&A, 336, 44 

Poutanen J., Vurm I., 2010, ApJS, 189, 286 

Pratt G. W., Arnaud M., 2002, A&A, 394, 375 

Prokhorov D. A., Colafrancesco S., 2012, MNRAS, 424, L49 

Prokhorov D. A., Colafrancesco S., Akahori T., Million E. T., Nagataki S., 

Yoshikawa K, 201 1, MNRAS, 416, 302 
Reese E. D., Carlstrom J. E., Joy M., Mohr J. J., Grego L., Holzapfel W. L., 

2002, ApJ, 581,53 
Rephaeli Y, 1995a, ARA&A, 33, 541 
Rephaeli Y, 1995b, ApJ, 445, 33 
Rudd D. H., Nagai D., 2009, ApJL, 701, L16 
Rudd D. H., Zentner A. R., Kravtsov A. V, 2008, ApJ, 672, 19 
Sazonov S. Y, Sunyaev R. A., 1998, ApJ, 508, 1 
Shimon M., Rephaeli Y, 2004, New Astronomy, 9, 69 
Sunyaev R. A., Norman M. L., Bryan G. L., 2003, Astronomy Letters, 29, 

783 

Sunyaev R. A., Zeldovich I. B., 1980, MNRAS, 190, 413 

Vikhlinin A., Kravtsov A., Forman W., Jones C, Markevitch M., Murray 

S. S., Van Speybroeck L., 2006, ApJ, 640, 691 
Williamson R. et al., 2011, ApJ, 738, 139 
Wright E. L., 1979, ApJ, 232, 348 

Zeldovich Y. B., Sunyaev R. A., 1969, Astrophysics and Space Science, 
4, 301 

Zemcov M. et al., 2012, ApJ, 749, 1 14 
Zemcov M. et al., 2010, A&A, 518, L16 

© 0000 RAS, MNRAS 000, 000-000 



